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ABSTRACT 

First-principles  calculations  based  on  spin-polarized  density  functional  theory  (DFT) 
and  the  generalized  gradient  approximation  (GGA)  have  been  used  to  study  the  adsorption 
of  nitromethane  (NM)  and  1,1 -diamino-2, 2-dinitroethylene  (FOX-7)  molecules  on  the 
Al(lll)  surface.  The  calculations  employ  (3x3)  aluminum  slab  geometries  and  3D 
periodic  boundary  conditions.  Based  on  these  calculations,  we  have  determined  that  both 
dissociative  and  nondissociative  adsorption  mechanisms  are  possible,  depending  on  the 
molecular  orientation  and  the  particular  surface  sites  involved.  In  the  case  of  dissociative 
chemisorption,  O  abstraction  by  A1  surface  atoms  is  seen  to  be  the  dominant  mechanism. 
The  dissociated  oxygen  atom  forms  strong  Al-0  bonds  with  the  neighboring  A1  sites 
around  the  dissociation  sites.  Additionally,  the  radical  species  obtained  as  a  result  of 
oxygen  atom  elimination  remains  bonded  to  the  surface.  In  some  instances,  both  oxygen 
atoms  of  the  nitro  group  dissociate  and  oxidize  the  aluminum  surface.  Finally,  for  the 
case  of  nondissociative  adsorption,  various  N-O-Al  bridge-type  bonding  configurations 
can  be  formed.  Based  on  the  data  provided  from  these  studies,  it  can  be  concluded  that 
oxidation  of  the  aluminum  surface  readily  occurs,  either  by  partial  or  complete 
dissociation  of  the  oxygen  atoms  from  the  NO2  group. 
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I.  Introduction 

Powderized  aluminum  has  long  been  used  as  an  energetic  ingredient  in  rocket 
propellant  formulations,  comprising  approximately  15-20%  of  conventional  ammonium 
perchlorate  solid  propellant.1  Its  primary  roles  are  to  increase  the  combustion 
exothermicity  and  the  regression  rate  of  solid  propellant  grains  and  to  enhance  the  blasting 
effect  of  explosives.  Such  properties  are  significantly  affected  by  the  size  of  the  A1 
particles.  For  example,  in  the  case  of  A1  nanopowders,  significant  improvements  in  the 
performance  of  some  energetic  materials  over  the  common  micron-size  A1  powders  have 
been  reported."" 

An  important  performance  issue,  particularly  for  nanoscale  aluminum  powders,  is 
the  formation  of  an  aluminum  oxide  overcoat  prior  to  combustion  which  inhibits  efficient 
burning.  When  the  thickness  pf  the  oxide  overcoat  is  large,  its  presence  severely  reduces 
the  potential  advantages  of  using  such  high  surface  area-to-volume  ratio  ultrafine 
aluminum  particles.  Consequently,  identification  of  new  methods  for  passivation  of  A1 
particles  to  reduce  the  size  of  the  oxide  layer,  without  incurring  a  significant  decrease  in 
performance,  are  being  considered.  For  example,  it  has  been  suggested  to  coat  the 
aluminum  particles  with  an  energetic  material  such  as  hexahydro-1 ,3,5-trinitro-.i>-triazine 
(RDX).4 

Currently,  no  experimental  evidence  is  available  to  clarify  the  efficiency  of 
passivation  of  A1  nanoparticles  when  coated  with  various  energetic  materials.  Moreover, 
a  clear  mechanism  to  describe  the  interaction  of  energetic  materials  (nitro  compounds  in 
particular),  with  A1  surfaces  is  lacking.  In  order  to  clarify  some  of  the  fundamental  issues 
related  to  the  role  played  by  energetic  materials  when  deposited  on  an  aluminum  surface, 
the  present  work  focuses  on  the  atomic-level  description  of  the  interactions  between  the 
energetic  compounds  nitromethane  (NM)  and  l,l-diamino-2,2-dinitroethylene  (FOX-7) 
with  the  Al(lll)  surface.  Specifically,  first-principles  quantum  chemical  calculations  are 
used  to  determine  the  chemisorption  properties  of  the  nitro-containing  compounds  of 
interest  for  energetic  materials  applications  on  aluminum  surfaces. 
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II.  Computational  Method 

The  calculations  performed  in  this  study  were  done  using  the  Vienna  ab  initio 
simulation  package  (VASP).5'7  This  program  evaluates  the  total  energy  of  periodically 
repeating  geometries  based  on  density-functional  theory  and  the  pseudopotential 
approximation.  The  electron-ion  interaction  is  described  by  fully  non-local  optimized 
ultrasoft  pseudopotentials  similar  to  those  introduced  by  Vanderbilt.8,  Periodic 
boundary  conditions  are  used,  with  the  one  electron  pseudo-orbitals  expanded  over  a 
plane- wave  basis  set.  The  expansion  includes  all  plane  waves  whose  kinetic  energy, 

r|2k2  /  2m  <  Ecut  where  k  is  the  wave  vector,  m  the  electronic  mass  and  Ecut  is  the 
chosen  cutoff  energy.  In  this  study  a  cutoff  energy  of  395  eV  is  chosen  which  ensures  the 
convergence  with  respect  to  the  basis  set. 

Calculations  were  performed  using  spin-polarized  generalized  gradient 
approximation  (GGA)  density  functional  theory  using  PW91  exchange-correlation 
functional.10  The  sampling  of  the  Brillouin  zone  was  performed  using  a  Monkhorst-Pack 
scheme.11  We  have  also  used  the  Methfessel-Paxton  technique12  with  a  smearing  of  a=0.2 
eV  in  order  to  minimize  the  errors  in  the  Hellmann-Feynman  forces  due  to  the  entropic 
contribution  to  the  electronic  free  energy.5,6  All  energies  are  extrapolated  to  T=0  K. 

The  minimization  of  the  electronic  free  energy  was  performed  using  an  efficient 
iterative  matrix-diagonalization  routine  based  on  a  sequential  band-by-band  residuum 
minimization  method  (RMM),5,6  or  alternatively  based  on  preconditioned  band-by-band 
conjugate-gradient  (CG)  minimization.18  The  optimization  of  different  atomic 
configurations  was  performed  by  conjugate-gradient  minimization  of  the  total  energy. 

The  A1  surface  is  represented  by  a  slab  model  with  periodic  boundary  conditions 
applied  in  all  three  directions.  Particularly,  a  3x3  supercell  with  four  layers  containing  36 
A1  atoms  is  used  to  study  the  adsorption  of  various  molecular  systems.  The  slabs  were 
separated  by  10  A  and  15  A  of  vacuum  for  the  cases  of  NM  and  FOX-7  molecules, 
respectively.  In  calculations  of  molecular  adsorption  on  the  surface,  we  have  relaxed  all 
atomic  positions  of  the  molecule  as  well  as  the  A1  atoms  located  in  the  first  layer  of  the 
slab. 
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III.  Results  and  Discussion 

A.  Test  Calculations  for  Bulk  A1  and  isolated  NM  and  FOX-7  molecules 

A  number  of  tests  have  been  initially  performed  to  verify  the  accuracy  of  the 
method  when  applied  to  bulk  A1  and  to  the  isolated  NM  and  FOX-7  molecules,  and  to 
clarify  different  technical  aspects  such  as  the  optimum  cutoff  energy  for  calculations. 

For  bulk  Al,  we  have  tested  for  convergence  the  k-point  sampling  density  as  well 
as  the  kinetic  energy  cutoff.  In  these  calculations  a  Monkhorst-Pack  grid  of  12x12x12  has 
been  used,  leading  to  56  k-points  in  the  irreducible  Brillouin  zone.  In  order  to  deteimine 
the  equilibrium  bulk  parameters  of  Al  we  have  uniformly  scaled  the  lattice  vectors  and 
performed  energy  calculations  as  a  function  of  unit  cell  volume.  The  calculated  data  were 
then  fit  with  the  Murnaghan  equation  of  state.14  The  calculated  lattice  constant  was  found 
to  be  4.044  A,  which  differs  by  -0.14%  from  the  corresponding  experimental  value  of 
4.050  A.15  The  corresponding  bulk  modulus  and  its  pressure  derivative  were  found  to  be 
B o=72. 2  GPa  and  B0’=4.1 6, |; respectively.  These  values  differ,  respectively,  by  -4.87% 
and  -2.57%  from  the  corresponding  values  of  75.9  GPa  and  4.27,  estimated  in  Ref.  16 
based  on  the  reported  elastic  constants  of  Al.  These  values  indicate  that  the  present  set  of 
pseudopotentials  is  able  to  provide  a  very  good  representation  of  the  structural  properties 
of  bulk  Al. 

An  equally  good  representation  has  been  observed  for  the  geometric  parameters  of 
the  isolated  NM  and  FOX-7  molecules.  For  example,  based  on  optimizations  of  the 
isolated  NM  molecule  in  a  cubic  box  of  12x12x12  A3,  we  have  determined  the  following 
equilibrium  bond  lengths:  r(C-N)=1.497  A,  r(N-O)=1.240  A,  and  r(C-H)=1.090  A,  and 
bond  angles  0(N-C-H)=1O6.6°  and  0(N-O-N)=125.6°.  These  values  are  very  close  to 
those  we  have  reported  previously17  for  the  isolated  NM  molecule  based  on  hybrid 
density  functional  theory  calculations  at  the  B3LYP/6-31G**  theoretical  level. 
Additionally,  the  calculated  values  are  also  close  to  the  experimental  data  for  gas  phase 
nitromethane.18  In  this  case  the  largest  deviation  of  1.3%  is  seen  for  N-0  bonds  (rexp(N- 
0)=  1.224  A). 

For  the  case  of  the  FOX-7  molecule,  we  have  shown  in  a  previous  study19  that 
plane-wave  DFT  calculations  using  the  VASP  code  provide  an  accurate  description  of  the 
structural  properties  of  this  molecule  as  well  as  the  corresponding  crystal  structure.  The 
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corresponding  analysis  will  not  be  repeated  here.  The  good  agreement  between  our 
calculated  properties  of  bulk  A1  and  the  isolated  NM  and  FOX-7  molecules  with 
experiment  and/or  previous  theoretical  predictions  made  us  confident  to  proceed  to  the 
next  step,  i.e.,  investigation  of  molecular  adsorption  on  the  Al(l  1 1)  surface. 

B.  Nitromethane  Adsorption  on  Al(lll)  Surface. 

Bl.  Geometries  and  Energies 

The  adsorption  of  NM,  which  is  a  prototype  for  larger  energetic  molecules  such  as 
hexahydro- 1 ,3,5-trinitro-5-triazine  (RDX)  or  octahydro- 1,3,5, 7-tetranitro- 1,3, 5,7- 
tetrazocine  (HMX),  on  the  Al(lll)  surface  was  examined  in  detail  for  several  distinct 
orientations  of  NM  relative  to  the  metallic  surface.  For  the  case  of  non-dissociative 
adsorption  configurations,  the  corresponding  adsorption  energy  was  calculated  according 
to  the  expression 

Eads— Em0]ec  Es]ab  -  E(molec+slab) 

where  Emoiec  is  the  energy  of  the  isolated  NM  molecule  in  its  equilibrium  position,  Esiab  is 
the  total  energy  of  the  slab  and  E(m0]ec+siab)  is  the  total  energy  of  the  adsorbate/slab  system. 
A  positive  EadS  corresponds  to  a  stable  adsorbate/slab  system.  The  energy  of  the  isolated 
NM  molecule  for  this  calculation  is  that  determined  for  the  isolated  NM  molecule  in  a 
cubic  cell  of  length  12  A.  The  same  Brillouin-zone  sampling  has  been  used  to  calculate 
the  energies  of  the  bare  slab  and  of  the  molecule-slab  systems. 

Representative  adsorption  configurations  are  depicted  in  Figures  1  and  2.  For  each 
configuration,  three  views  are  represented.  The  first  one,  denoted  with  the  index  1, 
represents  the  initial  molecular  configuration  at  the  start  of  the  optimization  process.  The 
second  and  third  insets,  denoted  with  indices  2  and  3  ,  respectively,  represent  side  and  top 
views  of  the  final  adsorption  configuration  after  full  relaxation  of  atomic  positions. 

In  Figure  1  we  present  the  adsorption  configurations  of  NM  when  initially  the  C-N 
bond  is  oriented  perpendicular  to  the  surface,  with  the  nitro  group  pointing  down  toward 
the  surface.  Panels  al)-a3)  illustrate  the  results  obtained  in  which  adsorption  of  NM 
takes  place  at  an  fee  site.  In  this  case,  various  Al-0  bonds  are  formed  with  equilibrium 
values  between  1.897  and  2.098A.  The  final  (i.e.,  optimized)  configuration,  denoted  as 
NM(I),  represented  in  panels  a2)  and  a3)  illustrates  marked  distortions  of  the  nitro  group 
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due  to  the  strong  interactions  of  the  oxygen  and  nitrogen  atoms  with  the  aluminum 
surface.  As  a  result,  the  N-0  bonds  are  significantly  stretched  to  1.455  and  1.476  A, 
respectively,  relative  to  the  gas  phase  equilibrium  values  of  1.240  A.  Additionally,  the 
nitro  group  and  the  C  atom  are  no  longer  coplanar  as  in  the  gas  phase,  but  in  this  state  the 
molecule  adopts  a  chair  conformation.  The  corresponding  torsional  angle  x(O-O-N-C) 
changes  from  the  gas  phase  value  of  180°  to  -1 1 1.7°  for  this  chemisorbed  state.  The  large 
geometric  changes  that  result  upon  adsorption  are  also  reflected  in  a  high  adsorption 
energy.  Based  on  Eq.  (1)  and  the  diagram  level  provided  in  Figure  3,  the  predicted 
adsorption  energy  is  43.3  kcal/mol. 

Panels  bl)-b3)  of  Figure  1  illustrate  the  case  in  which  the  initial  configuration  of 
the  NM  molecule  interacts  with  an  aluminum  on-top  site.  As  can  be  seen  from  the  b2)  and 
b3)  images,  adsorption  at  this  site  leads  to  complete  dissociation  of  one  of  the  oxygen 
atoms  from  the  nitro  group.  The  dissociated  O  atom  adsorbs  at  a  nearby  fee  site  with 
simultaneous  bonding  to  three  neighboring  A1  atoms  in  the  first  layer.  The  corresponding 
Al-0  bonds  are  almost  identical,  having  values  between  1.85  and  1.86  A.  The  remaining 
NM  fragment,  denoted  Ntr(I),  also  chemisorbs  strongly  to  the  surface,  forming  bonds  to 
aluminum  via  both  the  nitrogen  atom  and  the  remaining  oxygen  atom.  The  corresponding 
Al-N  bond  length  is  1.849  A  while  the  Al-0  bond  lengths  are  1.896  and  1.904  A.  As  in  the 
case  of  the  adsorption  at  the  fee  site,  the  remaining  (undissociated)  N-0  bond  of  the 
nitroso  species  is  significantly  stretched,  with  a  length  of  1.504  A. 

In  order  to  determine  the  strength  of  the  adsorption  for  this  nitroso  complex,  we 
have  used  a  relation  similar  to  that  given  in  Eq.  (1)  in  which  Esiab  was  taken  as  the  energy 
of  the  slab  and  dissociated  oxygen  atom  complex  and  Emoiec  was  taken  to  be  the  energy  of 
the  isolated  nitrosomethane  species.  The  specific  values  of  the  energetic  levels  involved 
are  shown  in  Figure  3.  Relative  to  the  sum  of  energies  of  the  surface  with  an  adsorbed 
oxygen  atom  and  of  the  isolated  nitrosomethane  species  (see  EAi+o(a)+Ntr(g)  in  Figure  3a), 
the  adsorbed  nitroso  complex  at  the  fee  site  has  an  adsorption  energy  EAi+o(a)+Ntr(i)  of 
127.6-60.5=67.1  kcal/mol.  This  result  indicates  that  even  after  cleavage  of  the  first  N-0 
bond,  the  nitroso  complex  also  interacts  strongly  with  the  A1  surface. 

Finally,  we  present  in  the  panels  cl)-c3)  the  results  of  adsorption  at  an  hep  surface 
site.  In  this  case,  both  oxygen  atoms  of  the  nitro  group  dissociate,  leading  to  formation  of 
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a  methylnitride  chemisorbed  species  (denoted  MeN(I)).  The  dissociated  O  atoms  occupy 
fee  surface  sites  and  form  strong  bonds  to  the  A1  surface.  The  corresponding  Al-0  bonds 
have  values  of  1.847,  1.847  and  1.863  A.  The  longest  Al-0  bonds  correspond  to  those  A1 
atoms  which  are  also  involved  in  simultaneous  bonding  to  dissociated  O  atoms  and  the  N 
atom  of  the  methylnitride  chemisorbed  species  (see  panel  c3  in  Figure  1).  Likewise,  the 
nitrogen  atom  in  the  original  nitro  group  also  forms  bonds  to  the  metal  surface  with  Al-N 
bond  lengths  that  vary  from  1.926  to  1.937  A. 

The  strength  of  the  interaction  of  methylnitride  with  the  surface  E(Ai+20(a)+MeN(i))- 
E(Ai+20(a))-E(MeN(g))  has  been  determined  with  respect  to  the  sum  of  energies  of  the  surface 
with  two  adsorbed  oxygen  atoms  (EAi+20(a))  nnd  that  of  isolated  methylnitride  chemisorbed 
species  (EmcN®).  The  corresponding  energy  diagram  is  given  in  Figure  3a.  Based  on 
these  values  we  have  determined  the  adsorption  energy  to  be  116.1  keal/mol  for 
methylnitride  to  the  surface. 

In  addition  to  the  vertical  adsorption  configurations  described  above,  we  have  also 
tested  the  case  in  which  NM  is  initially  oriented  with  the  methyl  group  pointed  toward  the 
surface.  In  this  case,  chemisorption  does  not  occur  since  the  methyl  group  does  not 
interact  with  the  surface  aluminum  atoms.  Such  configurations  do  not  lead  to  stable 
adsorption  configurations  and  we  have  not  investigated  them  further. 

Beside  the  vertical  adsorption  configurations,  we  have  also  analyzed  the  case  in 
which  the  NM  molecule  is  initially  oriented  parallel  to  the  metallic  surface.  We  have 
examined  three  cases:  the  N  atom  of  NM  positioned  above  an  A1  surface  atom  (see  panel 
al.  Figure  2),  above  an  hep  site  (see  panel  bl,  Figure  2),  and  above  an  fee  site  (see  panel 
cl,  Figure  2).  Our  calculations  indicate  that  during  the  optimization  of  the  atomic 
positions  started  from  these  parallel  configurations,  the  NM  molecule  rotates  in  order  to 
maximize  the  interaction  of  the  nitro  group  with  the  aluminum  surface.  As  a  result,  the 
molecule  initially  tilts  such  that  the  oxygen  atoms  points  toward  the  surface.  Additionally, 
depending  on  the  particular  molecular  orientation,  adsorption  will  take  place  either 
molecularly  or  dissociatively. 

Panels  al)-a3)  of  Figure  2  illustrate  a  case  in  which  dissociation  does  not  take 
place,  denoted  as  NM(II).  In  this  case,  rather  large  molecular  deformations  take  place. 
As  illustrated  in  panel  a2)  of  Figure  2,  the  molecule  adopts  a  chair  configuration  with  the 


torsional  angle  t(0-0-N-C)=115.6°.  The  molecular  C-N  axis  is  tilted  by  28°  relative  to 
the  surface  normal.  As  a  result  of  intramolecular  deformation,  both  the  N  and  the  O  atoms 
of  the  nitro  group  participate  in  the  bonding  to  the  A1  surface.  In  this  case,  the  Al-0  bond 
lengths  are  1.840  A  while  the  Al-N  bond  length  is  2.111  A.  The  strong  interaction  with 
the  metallic  surface  is  also  reflected  by  the  stretched  N-0  bonds  of  the  NM  molecule, 
which  are  1.429  A,  0.19  A  larger  than  those  in  the  gas-phase  molecule.  For  the 
chemisorbed  configuration  with  the  N  atom  positioned  above  an  hep  site,  we  find  a 
binding  energy  ENm(ii)  of  53.8  keal/mol,  which  is  slightly  larger  than  that  found  for  the 
undissociated  state  in  which  the  optimization  was  started  from  a  vertical  configuration  of 
NM  (Enmcii))- 

In  the  case  where  the  initial  configuration  corresponds  to  the  N  atom  above  an  hep 
site,  adsorption  on  the  A1  surface  leads  to  dissociation  of  one  of  the  N-0  bonds  (see  panels 
b2  and  b3  of  Figure  2),  with  formation  of  chemisorbed  nitrosomethane,  denoted  as 
Ntr(II).  The  dissociated  0  atom  is  adsorbed  at  an  fee  site  and  forms  Al-0  bonds  with 
lengths  in  the  range  1.821-1.854  A.  The  remaining  O  atom  of  the  nitro  group  bonds  to 
two  surface  A1  atoms  (r(Al-0)=  1.877  A  and  1.938  A),  while  the  N  atom  binds  at  a 
distance  of  1.865  A  to  an  A1  surface  atom.  As  before,  an  elongated  N-0  bond  (1.523  A)  is 
observed.  The  corresponding  binding  energy  of  this  nitrosomethane  chemisorbed  species 
(ENtr(ii))  is  found  to  be  60.1  keal/mol,  close  to  the  value  obtained  for  Ntr(I). 

Finally,  in  the  case  where  the  N  atom  is  initially  positioned  above  an  fee  site 
(panel  cl,  Figure  2),  adsorption  leads  to  dissociation  of  both  O  atoms.  As  a  result,  the  O 
atoms  form  strong  Al-0  bonds  (see  panels  c2)  and  c3)  of  Figure  1)  while  the  remaining 
methylnitride  chemisorbed  species,  denoted  MeN(II),  adsorbs  perpendicular  to  the  surface 
(see  panel  c2)  of  Figure  1).  It  is  important  to  note  that  the  methylnitride  chemisorbed 
species  also  binds  strongly  to  the  surface.  Based  on  the  energy  values  represented  in 
Figure  3  we  have  obtained  an  adsorption  energy  of  116.2  keal/mol  for  the  methylnitride 
chemisorbed  species.  This  radical  is  bonded  to  the  surface  with  Al-N  bond  lengths  of 
1.924,  1.936  and  1.936  A. 

Overall,  these  results  indicate  that  NM  can  undergo  both  molecular  and 
dissociative  adsorption  when  deposited  on  the  A1  (111)  surface.  The  fact  that  dissociation 
of  the  nitro  group  on  the  A1  surface  is  observed  in  our  simple  energy  minimizations 
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suggests  that  the  activation  energies  of  these  processes  are  minimal.  Moreover,  we  note 
that  in  addition  to  the  formation  of  strong  Al-0  bonds,  the  dissociation  products 
nitrosomethane  or  methylnitride  also  strongly  bind  to  the  surface  through  interactions  of 
either  N  or  O  atoms. 


B2.  Electron  Localization  Function  and  Charge  Distribution 

A  useful  tool  to  investigate  the  nature  of  the  chemical  interaction  of  NM  with  the 
Al(l  1 1)  surface  is  the  electron  localization  function  (ELF),  which  is  defined  as  21,22 


ELF(r) 


1 


r  D(r)  V 


(2) 


1  + 


vDh(r)y 

where  D(r)  is  the  Pauli  excess  kinetic  energy  density,  defined  as  the  difference  between 
the  kinetic  energy  density  and  the  von  Weizsacker  kinetic  energy  functional 


D(r)=jLvtv.p(P,iV^N^- 
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and  Dh(r)  is  the  kinetic  energy  of  the  homogeneous  electron  gas  for  a  density  equal  to  the 
local  density 


h  52mv  ' 

In  the  case  of  the  free  electron  gas  (Fermi  gas),  the  value  of  the  ELF  function  is  by 
definition  0.5.  Local  maxima  of  the  ELF,  also  called  localization  attractors,  correspond  to 
regions  where  electrons  form  pairs  of  opposite  spin.  Correspondingly,  in  the  case  of  fully 
localized  electrons  such  as  lone  electrons  or  electrons  paired  in  a  covalent  bond  with 
different  spin,  the  ELF  is  close  to  its  upper  bound  of  unity. 

Isosurfaces  of  ELF  can  be  calculated  and  used  to  help  explain  the  bonding 
structure  and  charge  localization  in  a  chemical  system.  A  plot  of  the  isosurface  for  a 
constant  ELF=0.7  value  in  the  case  of  isolated  NM  molecule  is  presented  in  Figure  4a. 
For  this  isosurface,  several  attractors  can  be  observed  in  the  figure:  three  corresponding  to 
C-H  bonds,  one  to  the  C-N  bond,  and  two  pairs  of  attractors  on  each  oxygen  atom 
corresponding  to  the  lone  pair  electrons.  We  notice  that  localization  is  strongest  in  the  C- 
H  bonds  and  in  the  lone  electron  pairs.  Figure  4  illustrates  the  ELF  isosurfaces 


10 


corresponding  to  adsorption  configurations  of  NM  presented  in  Figure  1,  namely  the 
nondissociated  state  NM(I)  (Figure  4b),  the  state  with  one  O  atom  dissociated  Ntr(I) 
(Figure  4c),  and  finally  the  state  with  both  O  atoms  dissociated  MeN(I)  (Figure  4d). 

For  the  nondissociative  adsorption  configuration  (Figure  4b),  we  observe  that 
kidney-shaped  distributions  on  O  atoms  change  to  ring  shapes,  corresponding  to  formation 
of  bonds  to  A1  atoms.  Also,  a  new  attractor  is  present  on  the  N  atom,  indicating  new 
bonds  to  the  substrate.  In  the  case  of  the  configuration  with  one  dissociated  O  atom 
(Figure  4c),  major  changes  are  seen  for  the  dissociated  O  and  N  atoms.  In  this  case, 
there  is  a  small  charge  localization  on  the  dissociated  0  atom.  This  indicates  an  ionic  type 
of  bonding.  Additionally,  the  charge  localization  on  the  N  atom  increases  relative  to  the 
non-dissociated  configuration,  as  indicated  by  the  larger  kidney  shape  attracted  around 
this  atom.  Finally,  in  the  case  when  both  O  atoms  of  NM  dissociate,  three  different 
attractors  are  present  around  the  N  atom  (see  Figure  4d).  These  attractors  have  lobes 
directed  towards  the  neighboring  A1  atoms.  This  corresponds  to  formation  of  new  bonds 
between  the  methylnitride  and  the  surface. 

A  common  tool  to  provide  a  semi-quantitiative  measure  of  charge  transfer  is  the 
Mulliken  population  analysis23  method  in  which  the  electronic  charge  is  partitioned 
among  the  individual  atoms.  This  partitioning  has  been  done  by  using  the  formalism 
developed  by  Segall  et  al.24  in  the  context  of  plane  wave  calculations  as  implemented  in 
the  CASTEP25  package.  Selective  variations  of  the  individual  atomic  Mulliken  charges 
on  O  and  N  atoms  of  NM  are  presented  in  Figure  5a  for  the  case  of  the  isolated  molecule 
(configuration  1)  as  well  as  for  the  three  adsorbed  states  depicted  in  Figure  1 
(configurations  2-4).  The  results  of  these  calculations  indicate  that  significant  charge 
transfer  takes  place  as  a  result  of  adsorption.  In  the  case  of  nondissociative  adsorption, 
there  is  a  significant  increase  in  the  negative  charge  of  the  O  atoms  of  about  -0.3e  relative 
to  gas  phase  NM.  A  corresponding  decrease  of  about  0.4e  takes  place  on  the  positively 
charged  N  atom.  Corresponding  to  this  transfer  of  electrons  to  the  NM  molecule,  there  is 
an  increase  in  the  positive  charge  of  the  Al  atoms  bonded  to  the  NM  molecule,  with  values 
in  the  range  0.38-0.5 le. 

The  process  of  charge  transfer  also  occurs  for  the  dissociative  adsorption 
configurations.  For  example,  in  the  case  in  which  one  (see  configuration  3  in  Figure  5a) 
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or  both  O  atoms  dissociate  (see  configuration  4  in  Figure  5a),  the  Mulliken  charges  of 
these  atoms  decrease  to  a  value  of  about  -le.  This  is  accompanied  by  a  corresponding 
decrease  of  the  charge  on  the  N  atom  from  0.5e  for  the  nitroso  complex  to  -l.Oe  for  the 
methylnitride  case.  Simultaneously,  there  is  a  significant  increase  in  the  positive  charge 
on  the  A1  atoms  involved  in  the  bonding  with  either  O  atoms  or  with  Ntr  or  MeN  radicals, 
with  values  ranging  from  0.5  to  0.88e.  For  the  A1  atoms  simultaneously  bonded  with  two 
O  atoms  or  with  an  O  and  a  N  atom,  some  larger  charges  in  the  range  0.80-0.88e  are 
observed. 

Overall,  these  results  indicate  that  there  is  a  significant  charge  transfer  from  the 
metallic  surface  to  the  NM  molecule  and  its  radicals.  Moreover,  the  amount  of  charge 
transferred  increases  with  molecular  dissociation.  These  results  are  consistent  with  the 
increase  of  the  binding  energy  for  molecular  NM  and  its  decomposition  fragments  to  the 
surface. 

C.  FOX-7  Adsorption  on  Al(lll)  Surface. 

1 1 

Cl.  Geometries  and  Energies 

The  strong  interactions  of  the  nitro  group  of  NM  with  the  A1  surface  are  likely  to 
be  present  for  other  molecular  systems  containing  this  chemical  group.  In  order  to  test 
this  notion,  we  have  analyzed  the  adsorption  of  another  energetic  molecule,  namely  FOX- 
7.  This  organic  compound,  whose  synthesis  was  reported  relatively  recently,  is 
characterized  as  a  promising  high  energy  density  material  with  superior  shock-sensitivity 
properties.  As  in  the  case  of  NM,  we  have  determined  representative  adsorption 
configurations  of  FOX-7  on  the  Al(lll)  surface  for  both  vertical  and  parallel 
configurations,  as  illustrated  in  Figures  6  and  7,  respectively.  In  each  case  we  present 
both  the  initial  molecular  configuration  at  the  start  of  the  optimization  process  as  well  as 
the  final  state  obtained  by  full  optimization  of  the  atomic  coordinates. 

Figure  6  shows  three  such  sets  of  configurations  of  FOX-7  interacting  with  the 
Al(lll)  surface.  In  panel  al)  the  approximately  planar  FOX-7  molecule  is  initially 
positioned  with  the  C=C  axis  perpendicular  to  the  surface.  One  of  the  O  atoms  of  the  nitro 
groups  is  positioned  approximately  on  top  of  an  Al  atom  while  the  C  atoms  are  positioned 
above  an  fee  surface  site.  An  additional  perpendicular  configuration,  with  the  C  atoms  of 
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the  FOX-7  molecule  positioned  above  an  hep  site,  is  shown  in  panel  bl).  Finally,  in  panel 
cl)  the  C=C  molecular  axis  is  tilted  relative  to  the  surface  normal  such  that  both  O  atoms 
of  the  same  nitro  group  interact  with  neighboring  A1  atoms. 

From  the  results  presented  in  panels  a2)  and  a3)  of  Figure  6,  we  observe  that 
adsorption  takes  place  molecularly  with  binding  of  one  of  the  0  atoms  of  each  nitro  group 
to  an  A1  atom  underneath.  We  denote  this  configuration  as  FOX7(I).  The  Al-0  distances 
are  1.868  A.  For  this  configuration,  the  major  changes  of  the  internal  bonds  are  seen  for 

o 

the  N-0  bonds  involved  in  bonding  to  the  surface,  which  are  stretched  to  about  1.35  A. 
Additionally,  the  corresponding  A1  atoms  where  bonding  takes  place  are  raised  upwards 
relative  to  the  other  surface  atoms  by  about  0.33  A.  For  this  chemisorbed  state  (denoted 
as  Eai+fox7(I)  in  Figure  6b),  the  predicted  binding  energy  is  28.4  kcal.mol.  This  value  is 
smaller  than  the  one  obtained  for  the  undissociated  molecular  configuration  of  NM(I)  (see 
panel  a2),  Figure  1).  This  difference  can  be  understood  based  on  the  number  of  direct 
bonds  formed  by  the  molecule  to  the  surface.  Indeed,  we  have  noted  that  in  the  case  of 
NM  there  is  a  significant  deformation  of  the  nitro  group,  which  adopts  a  chair-type 
conformation  such  that  both  ,0,  and  N  atoms  of  the  nitro  group  are  involved  in  the  bonding 
to  the  surface.  For  the  FOX-7  molecule,  however,  only  the  O  atoms  participate  in  direct 
bonding,  while  the  rest  of  molecular  structure  does  not  undergo  any  significant 
deformation.  The  lack  of  deformation  of  the  FOX-7  molecule  may  be  due  in  part  to  the 
intramolecular  N-H...0  hydrogen  bonds  between  the  amino  and  nitro  groups.  These 
intramolecular  interactions  tend  to  maintain  the  skeletal  structure  of  the  molecule. 

In  the  case  where  adsorption  is  initiated  in  a  vertical  configuration  with  the  C 
atom  above  an  hep  site  (see  panel  bl  in  Figure  6),  there  is  dissociative  chemisorption  of 
the  FOX-7  molecule  (see  panels  b2)  and  b3),  Figure  6),  leading  to  oxygen  abstraction  by 
the  surface  A1  atoms.  The  dissociated  0  atom  binds  to  three  neighboring  A1  atoms,  with 
Al-0  bond  lengths  in  the  range  1.833  tol.858  A.  The  remaining  nitroso-FOX-7  complex 
(denoted  as  NFOX7(I))  tilts  by  34°  relative  to  the  surface  normal,  thus  allowing  the 
interaction  of  both  nitro  and  nitroso  groups  with  the  surface.  In  this  case,  as  can  be  seen 
in  panel  b3)  of  Figure  6,  both  the  N  and  O  atoms  of  the  nitroso  group  interact 
simultaneously  with  one  of  the  surface  atoms,  denoted  All,  forming  bonds  of  lengths  r(N- 
Ah)=1.906  A  and  r(0-Ah)=1.792  A.  As  a  result  of  the  strong  interaction  with  the 
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molecule,  this  surface  Ah  atom  is  pulled  up  by  1.17  A  relative  to  the  surface  plane.  A 
similar  upward  displacement  is  seen  for  the  other  A1  atom,  denoted  Ah,  which  is  bonded 

o 

to  the  undissociated  nitro  group.  This  atom  is  moved  vertically  by  0.81  A. 

A  more  quantitative  measure  of  the  interaction  between  the  dissociated  FOX-7 
fragment  and  the  Al(lll)  surface  has  been  obtained  by  evaluation  of  the  adsorption 
energy.  Relative  to  the  energy  of  the  surface  with  an  adsorbed  oxygen  atom  and  that  of 
the  isolated  nitroso-FOX-7  radical  (EAi+o(a)+NFox7(g)),  the  adsorption  energy  is  69.1 
kcal/mol.  The  calculated  binding  energy  of  this  radical  is  similar  to  the  one  of  the  nitroso 
complex  of  67.1  kcal/mol  in  the  case  of  nitromethane  adsorption.  As  in  the  NM  case,  we 
find  that  after  cleavage  of  the  first  N-0  bond,  the  nitroso  complex  interacts  strongly  with 
the  A1  surface. 

Finally,  in  the  case  when  the  initial  configuration  of  the  FOX-7  molecule  is  in  a 
vertical  plane  but  with  the  C=C  axis  tilted  relative  to  the  surface  normal  as  illustrated  in 
panel  cl  of  Figure  6,  we  observe  that  adsorption  on  the  surface  leads  to  dissociation  of 
both  O  atoms  of  the  same  nitro  group,  with  formation  of  strong  Al-0  bonds.  The 
remaining  FOX-7  fragment  is  bonded  to  the  surface  through  the  N-atom  end  (r(N- 
Al)=1.878  A  and  1.916  A)  and  one  O  atom  of  the  other  nitro  group  (r(0-Al)=1.853  A). 
The  A1  atom  (denoted  as  All  in  panel  c3)  in  Figure  6)  simultaneously  bound  to  both  N 
and  O  atoms  of  the  FOX-7  fragment  is  displaced  upwards  relative  to  the  surface  plane  by 
1.38  A.  For  this  FOX-7  radical  (denoted  MFOX7(I)),  we  have  estimated  an  adsorption 
energy  of  160  kcal/mol.  This  binding  energy  is  higher  than  the  one  of  116.1  kcal/mol 
obtained  in  the  case  of  methylnitride  chemisorbed  species.  The  increase  in  binding  can  be 
understood  based  on  the  increase  in  the  number  of  direct  bonds  formed  between  the  FOX- 
7  radical  and  the  surface.  Indeed,  as  illustrated  in  Figure  6,  panel  c2,  the  FOX-7  radical  is 
bonded  to  the  surface  through  both  the  N  atom  of  the  dissociated  nitro  group  as  well  as 
through  one  of  the  oxygen  atoms  of  the  other  undissociated  nitro  group.  In  the 
methylnitride  case,  the  bonding  was  only  through  the  N  atom. 

In  Figure  7  we  present  the  results  of  atomic  optimizations  when  the  FOX-7 
molecule  is  initially  placed  in  a  plane  parallel  to  the  surface.  Two  configurations  have 
been  considered:  first,  (see  panel  al,  Figure  7)  one  O  atom  of  a  nitro  group  interacts 
initially  with  a  single  Al  atom  of  the  surface;  second,  (see  panel  bl,  Figure  7)  an  O  atom 
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of  a  nitro  group  interacts  simultaneously  with  two  neighboring  A1  surface  atoms.  As  in 
the  case  of  NM  adsorption  started  from  a  parallel  configuration,  the  FOX-7  molecule  first 
rotates  to  maximize  the  interaction  of  the  nitro  groups  to  the  surface.  Then,  depending  on 
the  particular  molecular  orientation,  adsorption  takes  place  either  molecularly  (see  panels 
a2  and  a3,  Figure  7)  or  dissociatively  (see  panels  b2  and  b3,  Figure  7). 

For  the  first  case,  shown  in  panels  a2  and  a3  of  Figure  7  (denoted  FOX7(II)),  the 
C=C  bond  is  tilted  by  about  60°  relative  to  the  normal  to  the  surface  and  all  four  0  atoms 
bind  independently  to  the  surface,  with  Al-0  bond  lengths  in  the  range  1.83-1.84  A. 
Additionally,  the  C  atoms  bonded  to  the  two  nitro  groups  forms  a  puckered  structure  (see 
panel  a3,  Figure  7)  with  torsional  angles  t(O-O-N-C)  in  the  range  118-122°.  The  values 
of  the  same  angles  for  the  gas  phase  structure  are  around  178°.  The  adsorption  energy  of 
this  configuration  is  69.5  kcal/mol.  The  increase  in  the  binding  energy  relative  to  the 
vertical  configuration  described  earlier  (see  panel  a2,  Figure  6)  is  clearly  due  to  the  fact 
that  now  there  are  four  Al-0  bonds  instead  of  two  bonds  as  seen  previously. 

Finally,  adsorption  initiated  from  a  parallel  configuration  with  simultaneous 
bonding  of  one  of  the  O  atoms  to  several  A1  surface  atoms,  results  in  dissociation  of  the 
corresponding  N-0  bond  (see  panels  b2  and  b3  of  Figure  7).  The  FOX-7  radical,  denoted 
NFOX7(II)  again  tilts  relative  to  the  surface  plane  such  that  the  C=C  bond  makes  an 
angle  of  about  57°  with  the  surface  normal.  The  dissociated  O-atom  binds  to  three  A1 
surface  atoms  while  the  other  O  atom  (denoted  with  Oi  in  Figure  5,  panel  b2)  of  the  nitro 
group  is  significantly  pulled:  toward  the  H  atom  of  the  amino  group  (denoted  with  Hi  in 
the  same  figure).  As  a  result,  the  corresponding  intramolecular  hydrogen  bond  distance 
r(O-H)  decreases  significantly  from  the  gas  phase  value  of  1.738  A  to  1.263  A.  The 
binding  energy  of  this  FOX-7  radical  is  84.9  kcal/mol. 

As  in  the  case  of  NM,  adsorption  of  FOX-7  on  the  Al(lll)  surface  can  lead  to 
both  molecular  and  dissociative  configurations.  For  the  dissociative  cases,  the  O  atoms  of 
the  nitro  group  embed  into  the  surface  and  form  strong  Al-0  bonds.  Additionally,  FOX-7 
fragments  bond  strongly  to  the  surface  through  either  N  or  O  atoms. 

C2.  Electron  Localization  Function  and  Charge  Distribution 
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A  comparison  of  the  changes  in  the  ELF  isosurfaces  for  gas  phase  FOX-7  and  for 
various  adsorption  configurations  is  presented  in  Figure  8.  For  the  isolated  FOX-7 
molecule  (see  Figure  8a),  several  attractors  can  be  identified.  For  example,  the  attractors 
corresponding  to  C=C,  C-N  and  N-H  bonds  are  clearly  evident  in  this  figure.  In  addition, 
the  attractors  corresponding  to  the  lone  pairs  of  electrons  on  the  N  atoms  of  the  amino 
groups  as  well  as  on  the  O  atoms  of  the  nitro  groups  are  also  present.  In  Figure  8  we  also 
illustrate  in  panels  b)-d)  the  ELF  isosurfaces  corresponding  to  the  adsorption 
configurations  of  FOX-7  shown  in  Figure  6,  namely  the  molecular  state  (Figure  8b),  the 
state  with  one  O  dissociated  (Figure  8c)  and  the  state  with  both  O  atoms  dissociated 
(Figure  8d). 

As  in  the  case  of  molecular  adsorption  of  NM,  molecular  adsorption  of  FOX-7 
induces  significant  changes,  in  the  charge  distributions  on  the  O  atoms.  In  this  case,  the 
kidney-shaped  distribution  seen  in  the  gas  phase  on  O  atoms  changes  to  a  ring  shape 
corresponding  to  formation  of  bonds  to  A1  atoms.  In  the  case  of  the  configuration  with 
one  dissociated  O  atom  (Figure  8c),  the  most  pronounced  changes  in  charge  distribution 
are  seen  for  the  O  and  N  atoms  involved  in  dissociation.  There  is  little  localization  of 
charge  on  the  dissociated  O  atom,  indicative  of  ionic  bonding.  However,  charge 
localization  on  the  N  atom  of  the  nitro  group  where  dissociation  takes  place  increases 
significantly  such  that  now  a  new  attractor  is  present  on  this  atom  (see  Figure  8c). 
Finally,  in  the  case  where  both  O  atoms  of  FOX-7  dissociate,  two  different  attractors  are 
now  present  around  the  N  atom,  as  shown  in  Figure  8d.  These  lobes  are  directed  towards 
the  two  neighboring  A1  atoms  and  correspond  to  formation  of  new  bonds  between  the 
FOX-7  radical  and  the  metallic  surface. 

As  in  the  case  of  NM  adsorption,  the  most  important  charge  variations  take  place 
on  the  nitro  group  of  the  FOX-7  molecule  involved  in  direct  bonding  to  the  surface.  The 
corresponding  variation  of  the  Mulliken  charges  on  the  corresponding  the  N  and  O  atom 
are  presented  in  Figure  5b..  Here  we  have  considered  the  case  of  isolated  molecule 
(configuration  1)  as  well  as  of  the  chemisorbed  states  illustrated  in  Figure  6 
(configurations  2-4).  The  overall  trend  seen  for  the  FOX-7  molecule  is  similar  to  that 
observed  for  NM  (see  Figure  5a).  Namely,  as  a  result  of  either  molecular  or  dissociative 
adsorption  on  the  surface,  there  is  a  significant  electronic  charge  transfer  from  the  surface 
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to  the  molecule.  The  amount  of  charge  transferred  increases  in  the  order:  molecular 
adsorbed  state  (configuration  2)  -»  dissociative  state  with  one  dissociated  O  atom 
(configuration  3)  dissociative  state  with  both  O  atoms  dissociated  (configuration  4). 
As  shown  in  Figure  5b,  there  is  a  continuous  decrease  of  the  charges  on  both  N  and  O 
atoms  when  going  from  configuration  2  to  configuration  4.  In  the  case  of  configuration  4 
of  Figure  5b,  the  charges  on  dissociated  0  atoms  as  well  as  on  the  N  atom  of  the  nitro 
group  involved  in  dissociation  are  about  -le.  Corresponding  to  this  transfer  of  elections  to 
the  FOX-7  molecule,  there  is  a  corresponding  increase  in  the  positive  charge  of  the  A1 
atoms  bonded  either  to  the  dissociated  O  atoms  or  to  the  remaining  FOX-7  fragment.  For 
example,  the  A1  atoms  bonded  with  a  single  O  atom  have  charges  in  the  range  0.38-0.46e. 
This  positive  charge  increases  to  0.94e  on  the  A1  atoms  bonded  simultaneously  with  two 
O  atoms.  The  largest  charge  increase  is  seen  for  those  A1  atoms  that  are  bonded 
simultaneously  to  O  and  N  atoms,  namely  atoms  denoted  with  All  in  Figure  6,  panels  b2 
and  c2.  The  corresponding  charges  on  these  atoms  are  1.26e  and  1.41e.  Overall,  as  in  the 
case  of  NM,  the  increase  in  the  amount  of  charge  transferred  to  the  adsorbed  molecule  or 
its  dissociation  fragments  is  directly  correlated  to  the  increase  in  the  corresponding 
binding  energy. 

IV.  Conclusions 

The  interactions  of  NM  and  FOX-7  molecules  with  the  Al(lll)  surface  have 
been  investigated  based  on  optimizations  performed  using  spin-polarized  GGA  plane- 
wave  density  functional  theory  calculations.  The  main  objective  of  this  study  was  to 
determine  the  possible  reaction  mechanisms  at  the  molecule-metallic  surface  interface. 
The  results  of  our  calculations  indicate  a  common  propensity  for  oxidation  of  the 
aluminum  surface  by  the  oxygen-rich  nitro  groups  (NO2)  of  the  NM  and  FOX-7 
molecules.  Dissociation  of  the  nitro  groups  and  formation  of  strong  Al-0  bonds 
appears  to  be  a  common  mechanism  for  these  molecules.  Additionally,  we  have  found 
that  both  the  molecular  adsorption  as  well  as  the  dissociation  processes  takes  place 
with  significant  charge  redistribution.  These  large  charge  redistributions  have  been 
found  to  be  consistent  with  the  strong  binding  energies  for  the  dissociated  O  atoms  and 
for  the  dissociation  fragments  of  either  NM  or  FOX-7  molecules.  These  results  suggest 
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that  contact  between  nitro-containing  energetic  compounds  and  A1  metal  will  result  in 
rapid  oxidation  reactions. 
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FIGURE  CAPTIONS 


Figure  1.  Adsorption  configurations  of  NM  on  the  Al(lll)  surface  obtained  from  initial 
vertical  configurations:  (a)  N  above  a  fee  site  (NM(I));  (b)  N  above  on-top  site  (Ntr(I)); 
(c)  N  above  an  hep  site  (MeN(I)).  The  initial  configurations  are  depicted  in  panels  al-cl 
while  the  lateral  and  top  views  of  the  corresponding  final  configurations  are  shown  in 
panels  a2-c2  and  a3-c3,  respectively. 

Figure  2.  Adsorption  configurations  of  NM  on  the  Al(lll)  surface  obtained  from  initial 
parallel  configurations:  (a)  N  above  a  top  site  (NM(II));  (b)  N  above  an  hep  site  (Ntr(II)); 
(c)  N  above  a  fee  site  (MeN(II)).  The  initial  configurations  are  depicted  in  panels  al-cl 
while  the  lateral  and  top  views  of  the  corresponding  final  configurations  are  shown  in 
panels  a2-c2  and  a3-c3,  respectively. 

Figure  3.  Relative  energies  of  various  configurations  of  NM  (a)  and  FOX-7  (b)  adsorbed 
on  the  Al(l  1 1)  surface  as  described  in  text.  The  abbreviations  used  corresponding  to:  h  - 
horizontal  configuration,  v  -  vertical  configuration,  a  -  adsorbed  configuration,  g  -gas 
phase  configuration.  The  other  abbreviations  NM(I),  NM(II),  Ntr(I),  Ntr(II),  MeN(I), 
MeN(II),  FOX7(I),  FOX7(II),  NFOX7(I),  NFOX7(II),  MFOX7(I)  are  described  in  the 
text. 

Figure  4.  Comparison  of  the  isosurfaces  of  the  electron  localization  function  (ELF=0.7) 
for  (a)  NM  in  the  gas  phase;  (b)  NM  molecularly  adsorbed  (NM(I));  (c)  NM  dissociatively 
adsorbed  with  one  O  atom  dissociated  (Ntr(I));  (d)  NM  dissociatively  adsorbed  with  both 
O  atoms  dissociated  (MeN(I)).  Configurations  (b)-(d)  correspond  to  the  same  structures 
shownin  Figure  1. 

Figure  5.  Variation  of  the  Mulliken  charges  on  the  individual  atoms  of  the  nitro  group  for 
(a)  NM  and  (b)  FOX-7.  In  both  cases  configuration  1  corresponds  to  gas  phase  molecules, 
while  configurations  2-4  correspond  to  the  molecular  adsorption  state,  dissociative 
adsorption  state  with  one  O  dissociated,  and  dissociative  adsorption  state  with  both  0 
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dissociated,  respectively.  These  configurations  are  those  depicted  in  Figure  1  for  NM  and 
in  Figure  6  for  FOX-7. 

Figure  6.  Adsorption  configurations  of  FOX-7  on  the  Al(lll)  surface  obtained  from 
initial  vertical  configurations:  (a)  C  atom  above  a  fee  site  (FOX7(I));  (b)  C  above  an  hep 
site  (NFOX7(I));  (c)  N  above  an  hep  site  (MFOX7(I)).  The  initial  configurations  are 
depicted  in  panels  al-cl  while  the  lateral  and  top  views  of  the  corresponding  final 
configurations  are  presented  in  panels  a2-c2  and  a3-c3,  respectively. 

Figure  7.  Adsorption  configurations  of  FOX-7  on  the  Al(lll)  surface  obtained  from 
initial  parallel  configurations:  (a)  O  atom  interacts  with  a  single  A1  atom  (FOX7(II));  (b) 
O  atom  interacts  simultaneously  with  two  A1  atoms  (MFOX7(II)).  The  initial 
configurations  are  depicted  in  panels  al-cl  while  the  lateral  and  top  views  of  the 
corresponding  final  configurations  are  presented  in  panels  a2-c2  and  a3-c3,  respectively. 

Figure  8.  Comparison  of  the  isosurfaces  of  the  electron  localization  function  (ELF=0.7) 
for  (a)  FOX-7  in  the  gas  phase;  (b)  FOX-7  molecularly  adsorbed  (FOX7(I));  (c)  FOX-7 
dissociatively  adsorbed  with  one  O  atom  dissociated  (NFOX7(I));  (d)  FOX-7 
dissociatively  adsorbed  with  both  O  atoms  dissociated  (MFOX7(I)).  Configurations  (b)- 
(d)  correspond  to  the  same  structures  as  shown  in  Figure  4. 
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